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Abstract 

o. 

Phase separation of binary fluids quenched by contact with cold external walls is considered. 
Navier-Stokes, convection-diffusion, and energy equations are solved by lattice Boltzmann method 
§ \ coupled with finite-difference schemes. At high viscosity, different morphologies are observed by 

varying the thermal diffusivity. In the range of thermal diffusivities with domains growing parallel 

£ 

to the walls, temperature and phase separation fronts propagate towards the inner of the system 
with power-law behavior. At low viscosity hydrodynamics favors rounded shapes, and complex 



o 



patterns with different lengthscales appear. Off-symmetrical systems behave similarly but with 
more ordered configurations. 
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I. INTRODUCTION 



When in a multi-phase system initially in a mixed state the temperature is decreased 
to values corresponding to a coexisting region of the phase diagram, domains of ordered 
phases start to form and grow with time. The process is called phase separation and is 
relevant for a large variety of systems [lj. In most of the cases studied theoretically, the 
temperature or other control parameters are assumed not depending on time and space, but 
are instantaneously set to their final values for coexistence. This assumption, reasonable in 
many situations, typically gives rise to a self-similar growth behavior with a characteristic 
domain size following a time power-law j^J. However, there are cases where the dynamics of 
the control parameter needs to be considered since it can greatly affect the morphology of 
domains. In binary alloys, for example, slow cooling is used to produce optimal sequences of 
a>ter„ate bands of different materia.* Q. In po.ymenc mixtures the possibility of controlling 
the demixing morphology by appropriate thermal driving has been studied in Refs. [5|, |6[; 
modulated patterns have been observed when a mixture is periodically brought above and 



below the critical point 



Q. 



Other worth examples of complexpattern formation due to 
the dynamics of the control parameters occur in crystal growth |8j, immersion-precipitation 
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[ill. 



membranes |9|, or in electrolyte diffusion in gels 

In this paper we study binary fluids quenched by contact with cold walls at temperatures 
below the critical value. The behavior of binary fluids in sudden quenches at homogeneous 
temperature is quite known j^, [l^] • For symmetric composition, the typical interconnected 
pattern of spinodal decomposition is observed. In the system here considered, phase separa- 
tion is expected to start close to the walls and develop in the inner of the system following the 
temperature evolution. The dynamics of this process and the role of the velocity field have 
not been explored too much, in spite of their relevance for many of the systems mentioned 
above. 

Two-dimensional studies of diffusive binary systems with cold sharp fronts propagating 
at constant speed have shown the formation of structures aligned on a direction de pen ding 



on the speed 
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161 ] . These results are also supported by theoretical analysis 
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Lamellar-like structures have been also found in numerical studies of two-dimensional off- 
symmetrical binary systems with the temperature following a fixed diffusive law 17|. In a 



model with the temperature dynamically coupled to the concentration field, point-like cold 



2 



sources have been shown to give rise to ring structures of alternate phases 18j . On the other 



hand, more usual morphologies have been found in cases with fixed thermal gradient [19], 
while complex phenomena such as sequential phase-separation cascades have been observed 
when the control parameter is slowly homogeneously changed j^ . The effects of full coupling 
between all thermo-hydrodynamic variables have been not considered sofar. 

The paper is organized as follows. In the next section the theoretical model and the 
numerical method are illustrated. The dynamics of our system is described by mass, mo- 
mentum, and energy equations with thermodynamics based on a free-energy functional 
including gradient terms. In Section III the results of our simulations are shown. We will 
explore the control parameter space by varying the viscosity and the thermal diffusivity. 
This will allow to analyze the differences with respect to the behavior of binary fluids in 
instantaneous quenching. The presentation will be focused on few cases typical for each 
regime. A final discussion will follow in Section IV. 

II. THE MODEL 

We consider a binary mixture with dynamical variables T,v, n, if which are, respectively, 
the temperature, the velocity, the total density, and the order parameter field being the 
concentration difference. Equilibrium properties are encoded in the free-energy 

(ij>(n,<p,T) + ±M\V<p\ 2 )dT (1) 



where 



^(n, <p,T) = e- k B T[n ln(n) - H—^) ~ ln (^^)] ( 2 ) 



with e = nksT +^(1 — ^2) being the bulk internal energy and the term in square brackets the 
mixing entropy. The gradient term in Eq. (1) is a combination of an internal ener gy g radient 
contribution proportional to K and of an entropic term proportional to — C [21| . hence 
M = K+CT. The system has a critical transition at /csT c = A/2 and the order parameter in 



the separated phases takes the values tp±(T) = ±y 3n 2 (T c /T — 1). The dynamical equations 
are given by 22] 

d t n = -d a (nv a ), (3) 
dtip = -d a {ipv a ) - 2d a J d w (4) 
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d t {nvp) = -d a (nv a vp) - d a (Jl a/3 - a a/3 ), 



(5) 



d t e = -d a {ev a ) - (II Q/3 - a a p)d a vp - d a J^, (6) 

where 3 d and 3 q are the diffusion and heat currents, H a p is the reversible stress tensor, 
o"a/3 = v(9 a V/3 + dpv a ) + (£ — 2r}/d)8 a pd y v^ is the dissipative stress tensor with (, rj being the 
bulk and shear viscosities, respectively, d the space dimension, and e = e + y|Vy?| 2 the total 
internal energy density also including gradient contributions. We have recently established 
the expressions for the pressure tensor Tl a p and chemical potential \i [23] following the 
approach of Ref . 2jJ . One finds 



n a/3 = (p- M^VV - M\ Vy\ 2 /2 - 2VvV • V(M/T)) 8 af} + Md a <pdp(p 



(7) 



where p = —ip + ndip/dn + ipdip/dip and \i = dip/d^T — TV • [(M/T)V</?]. Finally, in order 
to completely set up the dynamical system, phenomenological expressions for the currents 
are needed. As usually, one takes J d = -£nV(/i/T) + £i 2 V(l/T), J q = — £ 2 iV(/i/T) + 
£ 2 2V(1/T) where C a p is the positively defined matrix of kinetic coefficients with £ n = TV 
and £ 2 2 = T 2 k, T and k being the mobility and thermal diffusivity, respectively, assumed 



constant 



22|]. 



In order to solve Eqs. (EMJ) in d = 2 we have developed a hybrid lattice Boltzmann 
method (LBM) 24J427| where LBM [28( is used to simulate the continuity and Navier- 
Stokes equations ([3]) and (jSJ) while finite-difference methods are implemented to solve the 
convection-diffusion and the energy equations (jl]) and (J6]). LBM has been widely used to 



study mu 



ti-phase/component fluids 29] and, in particular, hydrodynamic effects in phase 
ordering 3j]]. It is defined in terms of a set of distribution functions, fi(r,t) with i = 
0, 1, 8, located in each site r at each time t of a D2Q9 (2 space dimensions and 9 lattice 
velocities) lattice where sites are connected to first and second neighbors by lattice velocity 
vectors of modulus |e$| = c (i = 1, 4) and |e*| = \/2c (i = 5, 8), respectively. The zero 
velocity vector e = is also included. The lattice speed is c = Ax/ At where Ax and At 
are the lattice and time steps, respectively. The distribution functions evolve according to 
a single relaxation time Boltzmann equation [31 ] supplemented by a forcing term {32! 



/<(r + e t At, t + At)- for, t) = - — [fi(r, t) - if (r, *)] + AtF t (r, t), 

T 



(8) 
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where r is the relaxation parameter, /? 9 are the equilibrium distribution functions, and Fj 
are the forcing terms to be properly determined. 

The total density and the fluid momentum are given by the following relations 

n = X;/i, nv = X;/iei + ^FAt, (9) 

where F is the force density acting on the fluid. The f^ q are expressed as a standard second 
order expansion in the fluid velocity v of the Maxwell-Boltzmann distribution functions 
331 ] . The forcing terms Fj in Eq. (jHJ) are expressed as a second order expansion in the lattice 
velocity vectors [34j. The continuity and the Navier-Stokes equations (ED and (0) can be 
recovered by using a Chapman-Enskog expansion when the Fj are given by 



/•; = i i - — ) 



e» - v Gj ■ v 
— — + ~^ e ' 



•F (10) 



with the force density F having components 

F a = d a (ncl)-d p U a ^ (11) 

c s = c/\/3 being the speed of sound in the LBM, w = 4/9, Ui = 1/9 for i = 1, ...,4, and 
= 1/36 for i — 5, 8. We observe that in this formulation the pressure tensor is inserted 
as a body force in the lattice Boltzmann equations. From the Chapman-Enskog expansion 
it comes out that £ = rj with 

n = nc 2 s At (-1 - iV (12) 



.At 2, 

On the other hand, a two-step finite difference scheme is used for the equations (jl]) and 
) (details on the implementation of Eq. (j4j) in the case of an isothermal LBM can be 



found in Ref. (27j). At walls, no-slip boundary conditions are adopted for the LBM |35|], the 
temperature is set to fixed values T& at the bottom wall and T u at the up wall, respectively, 
and neutral wetting for the concentration is adopted. This latter condition corresponds to 
impose a • Vip\ wa ii s = and a ■ V(V 2 ip)\ wa u s = 0, where a is an inward normal unit vector 
to the walls. These conditions together ensure a • V n\ wa u s = so that the concentration 
gradient is parallel to the walls and there is no flux across the walls. We have found this 
algorithm stable in a wide range of temperatures, viscosities and thermal diffusivities. With 



respect to thermal LBM for non-ideal fluids 36j where lattice Boltzmann equations are used 



to simulate the full set of macroscopic dynamical equations, the present model allows to 
reduce the number of lattice velocities thus speeding up the code and reducing the required 



memory 



27]. 



III. RESULTS AND DISCUSSION 



In the following we will explore the parameter space keeping fixed the values of K — 
0.003, C = 0, k B T c = 0.005, T = 0.1, and C\ 2 = £21 = 0. We will use lattices of size ranging 
from 256 x 256 to 1024 x 1024. We have considered different values of r\ and k. Before 
focusing on the cases representative of the various regimes, we will list all the runs we did 
in terms of dimensionless numbers. 

Common numbers used in hydrodynamics are the Reynolds and Peclet numbers Re and 
Pe. They are defined as Re = vL/u, where v = rj/n is the kinematic viscosity, Pe m d = vL/D 
for mass diffusion, where D is the mass diffusion coefficient, and Pe t d = vL/k for thermal 
diffusion. L and v are a typical length and velocity of the system. In phase separation 
L can be identified with the av erag e size of domains so that Re and Pe would depend on 
time (for a discussion see Ref. 37|). It is therefore more convenient for our purposes to 
introduce the Schmidt and Prandtl numbers Sc and Pr defined as Sc = vjD and Pr = v/k, 
where D = \a\T with a = (kfiT r /n)(T/T r — 1) being the coefficient of the linear term in the 
chemical potential \x (5), U^. Here T can be chosen as the value of the temperature at the 
walls. Table I contains a list of the runs we did, reported in terms of Sc and Pr. It is also 
useful to evaluate the Mach number Ma = \v\ max /c s where |v| ma:r is the maximum value of 
the fluid velocity during evolution. In all our simulations Ma is always much smaller than 
0.1 (see in the following), and the fluid results practically incompressible, as checked, with 
n ~ 1. For this reason we do not present in the paper any result about the time evolution 
of the total density n. 

First, as a benchmark for our method, we consider the relaxation of a single interface 
profile with k = 10~ 2 and rj = 0.167 (r = 1). This corresponds to a low viscosity regime as 
discussed in the following. We started the simulation with a sharp concentration step with 
values ip^(T b ) and <p + (T u ) and bulk temperature T/T c = 0.8 keeping fixed the temperatures 
Tt/T c = 0.8, T u /T c = 0.9 at the bottom and up walls (Fig. 1 (a)). The system reaches a 
stationary state with constant temperature gradient and concentration profile as in Fig. 1 (b). 
The numerical values of concentrations in the two bulk phases are in very good agreement 
with the analytical expression for ip±(T(r)) corresponding to the equilibrium values of T(r) 
shown in the related inset. This means that the concentration field p is in local equilibrium. 
The temperature of the up wall is then set to the same value of the temperature of the 
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bottom wall (Fig. 1 (c)). Then, as it can be seen in Fig. 1 (d), the system equilibrates 
at constant temperature with the expected concentration profile. Spurious velocities are of 
order 10~ 9 and result completely negligible. The test shows that stationary states are well 
reproduced by our algorithm. 



A. Diffusive regime 

We describe our results for phase separation. We first consider a case at very high viscosity 
with i] = 6.5 (r = 20) and symmetric composition (Runs 1-8). Here the effects of the velocity 
field are negligible. We set Tb/T c = T u /T c = 0.8 and initial bulk temperature above T c . As 
it can be seen in Fig. 2, for thermal diffusivities k > 10" 1 , usual isotropic phase separation 
is observed. In the range fc = 5x 10~ 4 4- 5 x 10 -2 , in spite of the neutral wetting condition 
on the boundaries, domains in the bulk have interfaces preferentially parallel to thermal 



fronts. For smaller values o 
agree with those of Refs. 13 



k domains grow perpendicularly to the walls. These results 



15 



16] in purely diffusive models where the same morphological 
sequence was found by decreasing the speed of cold fronts moving into a region with the 
mixed phase. However, also in absence of hydrodynamic effects, our case is different since 
the thermodynamics of the mixture is fully consistently treated and temperature fronts have 
no sharp imposed profile. 

We will now concentrate on cases at intermediate thermal diffusivities where domains are 
parallel to the walls and propagation fronts can be traced. Concentration and temperature 
configurations at successive times for k = 10~ 2 (Run 4a) are shown in Fig. 3 and Fig. 4, 
respectively. In this case it is Ma ~ 5 x 10 -5 . The temperature fronts have typical diffusive 
profiles which slowly relax to the equilibrium value imposed on the boundaries. In order 
to be quantitative, we defined yr(t) as the distance from the wall where the temperature 
assumes a fixed value (we chose T/T c = 0.88) and measured this quantity in simulations with 
large rectangular lattices. The solution of the diffusion equation with initial temperature 
T and fixed boundary value T w is (T(y,t) — T W )/(T — T w ) = erf[y/(2\/ki)] which implies 
yx/y/k ~ \ft. In the inset of Fig. 5 it is shown, in simulations with different fc, that yx 
follows the standard diffusion behavior. The time behavior of yx has been checked not 
depending on the specific value of the ratio T/T c in the range [0.8,1.0]; by considering a 
value of T such that T/T c < 1 allows to track the position of the temperature front for a 
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longer time interval. 

One can also consider the behavior of the fronts limiting the regions with separated phases, 
clearly observable in the first three snapshots of Fig. 3. Their position can be defined as the 
distance y v from the walls beyond which the condition V<£> ~ is verified everywhere. More 
precisely, we took y v as the point beyond which |V</?| < C with C = \/2 x 0.01; the value of 
C is chosen to match the maximum value of the fluctuations of \Vip\ in the initial disordered 
state, where \<p\ < 0.01. (In the last snapshot of Fig. 3 the two fronts propagating from up 
and down have come close each other and more usual phase separation occurs in the central 
region of the system.) We measured y v on rectangular lattices for different k and observed 
deviations from diffusive behavior (see Fig. 5). We found that y v grows by power law with 
an exponent depending on k. Our fits give y 9 ~ t 0,66 for k = 10~ 2 and exponents closer to 
1/2 for smaller k. We analyzed for different k possible variations of the typical values of 
fluid velocity but we did not find any. Therefore the change of the exponent of y v cannot be 
attributed to the velocity field. Even if y^ moves faster than yx and at long times it results 
Hip > Vt, we checked that the relation y v < yx c is always verified so that phase separation 
always occurs for T < T c . Since the phase separation is induced by the temperature change, 
one could have expected a similar behavior for y v and yx- The discrepancy could be related 
to the broad character of the temperature fronts which spreads the phase separated region. 
We also observed that the width of lamellar domains decreases at larger k, in agreement 



with Ref. [15 1 . 



B. Hydrodynamic regime 

At lower viscosities the evolution of morphology is very different in the range with inter- 
mediate values of thermal diffusivity. We will in particular illustrate in Fig. 6 the case with 
rj = 0.167 (r = 1) and k = 10~ 2 (Runs 19), for which we found Ma ~ 5 x 10~ 4 . This is 
the same thermal diffusivity of Fig. 3. At this viscosity hydrodynamics is relevant. Indeed, 
in instantaneous quenching at constant temperature and 77 = 0.167 we observed the domain 
growth exponent to assume the inertial value 2/3 (at odd with the diffusive high- viscosity 
value 1/3 ) [12]. The growth exponent was calculated by measuring the characteristic length 



defined by the inverse of the first momentum of the structure factor 



. The main effect 



due to hydrodynamics observable in Fig. 6 is that domains do not grow aligned with temper- 



S 



ature fronts as it occurs for the same thermal diffusivity at high viscosity. Circular patterns 



are stabilized by the flow 



12] and an example is given in Fig. 7. A similar picture occurs 



for other values of k here not reported (see Table I). On the other hand, the other thermal 
diffusivity regimes are less affected by hydrodynamics. When decreasing k, it is still possible 
to observe domains growing with interfaces normal to the walls as in the case at high viscos- 
ity (see Fig. 8 - Run 21b), while at larger k (Run 18) phase separation occurs isotropically 
like in an instantaneous quenching. 

The cases shown in Figs. 3 and 6 are typical of the high and low viscosity regimes. At 
intermediate values of 77 one can observe features common to the two above cases (see Fig. 9 
for 7] = 2.167 - Run 11a). Concerning the behavior of yr(t), we could not find relevant 
differences by varying r\ with respect to the case at high viscosity. 

Another effect induced by hydrodynamics is the formation of structures in the inner 
part of the system at earlier times than in the case at high viscosity (compare Fig. 3 and 
Fig. 6). In the inner region we can observe the typical interconnected pattern of spinodal 
decomposition but with a characteristic length-scale different from that of domains close 
to the walls. However, while the structures close to the walls are in local equilibrium, 
that is (p(r) = (p±(T(r)), in the middle of the system the concentration field is such that 
\tp\ < ip + (T(r)). A temporal regime characterized by the presence of domains with two scales 
was found in systems of different size (from 256 x 256 to 1024 x 1024) and k = 10~ 3 -v- 10~ 2 . 
In order to characterize the two scales we analyzed the behavior of the structure factor. In 
Fig. 10 the spherically averaged structure factor is shown at two consecutive times for a 
system having the same parameters of Fig. 6 and size L = 512. Two peaks are observable 
at each time that can be interpreted as related to the existence of two different length 
scales with one about twice longer than the other. The higher peak at smaller wave vector 
corresponds to the larger domains close to the walls while the other peak is related to the 
thinner domains in the inner of the system. At increasing times, the two peaks tend to 
merge. Due to this morphological evolution, in simulations at low viscosity, the position of 
the phase separation front y v could be measured only for a short time interval making not 
possible to determine the power-law behavior. 

Finally, we show results for systems with asymmetric composition. In Fig. 11 the evo- 
lution of two systems only differing for the value of viscosity is shown. Lamellar patterns 
prevail at high viscosity while circular droplets dominate at low viscosity (77 = 0.167). In the 
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latter case, again, two typical scales can be observed with thin tubes of materials connecting 
larger domains. The behavior of jjt is similar to that of the symmetric case. 

IV. CONCLUSIONS 

We have developed a numerical method for thermal binary fluids described by continu- 
ity, Navier-Stokes, convection-diffusion, and energy equations. We have studied quenching 
by contact with external walls, and we have shown how the pattern formation depends on 
thermal diffusivity, viscosity, and composition of the system. The evolution is very different 
from that observed in instantaneous homogeneous quenching. At high viscosity, different 
orientations of domains are possible. In an intermediate range of thermal diffusivities do- 
mains are parallel to the walls. The fronts limiting the regions with separated domains move 
towards the inner of the system with a power law behavior not always corresponding to that 
of the temperature fronts. At low viscosity, the velocity field favors more circular patterns, 
and domains are characterized by different length-scales close to the walls and in the inner 
of the system. Off-symmetrical mixtures give more ordered patterns. 

We conclude with two remarks on possible future directions of work. The first one 
concerns the Soret effect, which corresponds to have a mass diffusion current induced by 
thermal gradients. This effect can become relevant in quenching very close to the critical 
point where the ratio D T /D becomes large Here Dt is the thermal (mass) diffusion 
coefficient (Dt = C12/T 2 in our notation) and D is the mass diffusion coefficient defined 
at the beginning of Section III. In order to have a first idea on how the Soret effect can 
affect the pattern morphology, we considered a case with Dt/D = 20 corresponding to the 
highest values for this ratio reported in literature j^]. This would give Dt = 2 x 10 -3 , 
taking for D the value used in the runs of Section III. We run simulations for this case. We 
observed, in the intermediate range of thermal diffusivity and at high viscosity, the tendency 
of the system to exhibit more ordered lamellar patterns (parallel to the walls). At higher 
thermal diffusivity isotropic phase separation is found as usually, while at very low thermal 
diffusivity (k = 10~ 4 ), parallel patterns are found instead of perpendicular patterns. At 
low viscosity (we tested the case corresponding to that of Fig. 6) hydrodynamics continues 
to favor domains with more circular shape. We run also simulations with Dt = 10~ 4 , 
corresponding to a ratio Dt/D ~ 1, without finding relevant differences with the respect 
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to the case with Dt = 0. We also observe that the behavior of could depend on our 
choice for £12 and £21- A more comprehensive analysis of the Soret effect will be presented 
elsewhere. 

Finally, the morphology could be still richer in three dimensions, also due to the existence 
of more hydrodynamic regimes [2j, so that three-dimensional simulations would complete 
the picture given sofar. 
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Run 


Size 


Sc (xlO 3 ) 


Pr 


Symbol 


1 


512 


65 


12 


I 


2a, 2b 


512, 256 


65 


66 


I 


3 


512 


65 


129 


Pa 


4a, 4b 


512, 256 


65 


651 


Pa 


5a, 5b 


512, 256 


65 


1299 


Pa 


6 


256 


65 


6500 


Pa 


7 


512 


65 


65000 


Pe 


8 


256 


65 


650000 


Pe 


9 


256 


21.7 


22 


I 


10 


256 


21.7 


43 


I, Pa 


11a, lib 


512, 256 


21.7 


217 


I, Pa 


12 


256 


21.7 


2167 


Pe 


13 


256 


21.7 


21667 


Pe 


14 


256 


8.3 


8 


I 


15 


256 


8.3 


83 


I* 


16 


256 


8.3 


833 


I*, Pe 


17 


256 


8.3 


8333 


Pe 


18 


512 


1.7 


3 


I 


19a, 19b, 19c 


1024, 512, 256 


1.7 


17 


I* 


20a, 20b 


512, 256 


1.7 


167 


I* 


21a, 21b 


512, 128 


1.7 


1667 


Pe 



TABLE I: The first column indexes the simulation run, the second one is the linear size of the 
lattice, the third one is the Schmidt number (Sc), the fourth one is the Prandtl number (Pr). 
The last column is the symbol that identifies the kind of different observed patterns: I (isotropic 
morphology), Pa (domains parallel to the walls), Pe (domains perpendicular to the walls), I* 
(isotropic morphology with two lengthscales). The runs with two symbols exhibit patterns with 
common features to those corresponding to each symbol. 
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y y 

FIG. 1: Concentration and temperature (inset) profiles for an interface relaxation (see the text for 
explanation) . 



15 



FIG. 2: Typical configurations of the concentration field (p for symmetric composition at very high 
viscosity (77 = 6.5) with k = HT 1 , 1(T 2 , 1(T 5 (from left to right) at times t = 12.5 x 10 5 ;37.5 x 
10 5 ;300 x 10 5 , respectively, with lattice size 512 x 512, and T b /T c = T u /T c = 0.8. 
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FIG. 3: Configurations of concentration tp for composition 50/50 at times t = 7.5 x 10 5 ;22.5 x 
10 5 ; 37.5 x 10 5 ; 50 x 10 5 , at very high viscosity (77 = 6.5) with lattice size 512 x 512, T b /T c = T u /T c = 
0.8, and k = 10~ 2 . 
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FIG. 4: Configurations of the ratio T/T c for the same case and at same times of Fig. 3. Coordinates 
on the x and y axes are in lattice units and both of them are in the range [0, 512]. 
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io 6 t 10 7 

FIG. 5: Time behavior of y v at k = 10 -2 (empty symbols) and 1CT 3 (filled symbols) at very high 
viscosity with lattice size 128 x 2048. The straight line is a guide to the eye and has slope 2/3. 
Inset: Time behavior of yr/Vk at k = 10~ 2 (empty symbols) and 1CT 3 (filled symbols). The 
straight line has slope 1/2. 
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FIG. 7: Configuration of concentration ip for the case of Fig. 6 at time t = 11 x 10 5 with superim- 
posed the velocity field. 
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FIG. 8: Configuration of concentration ip at time t = 14 x 10 5 , low viscosity (77 = 0.167) as in 
Fig. 6, lattice size 128 x 128, T b /T c = T u /T c = 0.8, and k = 10~ 4 . 
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0.05 0.1 0.15 0.2 0.25 0.3 

k 

FIG. 10: Spherically averaged structure factor C(k) as a function of the wave vector modulus k 
for a system with the same parameters of Fig. 6 and size L = 512 at times t = 24 x 10 5 (empty 
symbols) and t = 39 x 10 5 (filled symbols), corresponding to the regime with two scales shown in 
Fig. 6. 
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FIG. 11: Configurations of concentration tp for composition 55/45 at times t = 4 x 10 5 ; 6 x 10 5 ; 11 x 
10 5 ; 16 x 10 5 , low viscosity (left column) and very high viscosity (right column) , lattice size 256x256, 
Tb/T c = T u /T c = 0.8, and k = 10~ 2 . Except for the composition, here the parameters are the same 
used in Figs. 6 and 3. 
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